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Abstract 

Backgrounds: Although height (H) has been considered the principal anthropometric variable governing lung function, the 
age-dependent differences in its influences on determining spirometric parameters (SPs) have not been conclusively 
investigated. Moreover, there has been no study centered on age-dependent effects of other anthropometric variables, 
including body weight (BW) and body fat mass (BFM) on SPs. In addition, the age-dependent influences of these 
anthropometric variables are anticipated to differ quantitatively between male and female participants. 

Methods: A total of 16,919 nonsmoking healthy Japanese adults (men: 6,116, women: 10,803) were partitioned into six 
groups stratified by gender and age at intervals of 20-years: young-, middle-, and advanced-age groups of either gender. 
Using a model in which a SP was described by a logarithmic additive function of age, H, BW, and BFM, we determined the 
partial regression coefficients of the respective anthropometric variables to predict the reference means of SPs, including 
FVC, FEV-,, FEV^FVC, PEF, FEF 50 , and FEF 75 , in the six groups. 

Results/Discussion: Although the impact of H on FVC and FEVt was relatively homogeneous irrespective of gender and age, 
its homogeneity faded for flow parameters, particularly in the female middle- and advanced-age groups, indicating that the 
age-dependent contribution of H to SPs was enhanced more in women. The impact of BW on SPs differed depending on 
age, and this effect was also more conspicuous for female participants. H and BW generally exerted positive effects on SPs, 
whereas BFM had negative effects. Opposite effects of BW and BFM were observed in the female middle-age group in 
particular. 

Conclusions: The effects of anthropometric variables on spirometric parameters are highly age-dependent, particularly in 
women, leading to the conclusion that the assumption of age-independent, constant partial regression coefficients of 
anthropometric variables while predicting the reference mean of a certain spirometric parameter may result in substantial 
errors. 
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Introduction 

Overwhelming numbers (more than 120 papers) of ethnic- 
specific regression equations for predicting reference means and/ 
or lower limits of normal (LLN) concerning various spirometric 
parameters in adults have been published over several decades [1]. 
Most of these equations were generated, however, by introducing 
age (A) and height (H) as explanatory variables. Incorporating the 
LMS (lambda, mu, and sigma) method with a smoothing function, 



Stanojevic et al. [2] recently demonstrated novel regression 
equations covering an entire range of age from preschool children 
to elderly persons in non-Hispanic white subjects. Extending the 
study of Stanojevic et al., the Global Lungs Initiative (GLI), an 
ERS Task Force, developed more global multi-ethnic regression 
equations with two explanatory variables, A and H, applicable 
across all ages from 3 to 95 years old for people of varied races 
[3-5]. The reason why most of the studies performed in this field 
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adopted only H as the explanatory anthropometric variable may 
be based on the fact that not only is the lung volume modeled most 
appropriately by the H (besides age), but, in addition, the goal of 
such studies is to develop a parsimonious equation that is easily 
handled in clinical situations [3]. Although this consideration is 
appreciated from a clinical point of view, it does not imply that 
other anthropometric variables play no role in deciding spiromet- 
ric parameters. In fact, numerous studies have been conducted to 
elucidate the effects of varied anthropometric variables other than 
H on spirometric parameters, though their impacts were not 
conclusively certified [6-19]. The common assumption in all the 
studies described above is that the effect of anthropometric 
variables, including H, body weight (BW), body mass index (BMI), 
or body fat mass (BFM) on a given spirometric parameter is 
independent of age. However, bodily shape or physique may 
change with aging. For instance, BW generally increases in middle 
age but decrease in advanced age [7] . H commonly decreases with 
increasing age. These facts indicate that the quantitative effect of a 
certain anthropometric variable on the decision of a spirometric 
parameter should not be taken to be constant but instead 
considered to change depending on age. Based on these 
considerations, we generated the following hypotheses in the 
present study: 1) spirometric parameters would be influenced by a 
variety of anthropometric variables, including H, BW, and BFM; 
2) their effects would differ depending on age, i.e., age-specific 
difference in contribution of an anthropometric variable to a 
spirometric parameter; and 3) the extent of age-dependent 
contribution of an anthropometric variable would differ qualita- 
tively and quantitatively between men and women, i.e., there are 
gender-specific differences in the effects of an anthropometric 
variable on a spirometric parameter. To test these hypotheses, we 
examined the relationship between anthropometric variables and 
spirometric parameters in six groups stratified by gender and age 
at intervals of 20 years. Each group consisted of thousands of 
healthy nonsmoking Japanese adults with a relatively narrow age 
span, i.e., 20 to 39 years old (young-age group), 40 to 59 years old 
(middle-age group) or more than 60 years old (advanced-age 
group). We believe that the investigation of groups with a narrow 
age range would be aid in extracting the age-specific contributions 
of anthropometric variables to spirometric parameter values. 

Methods 

Ethics Statement 

All participants provided written informed consent indicating 
they agreed their data could be used for clinical research. The 
participants were asked whether they agreed with the registration 
of their details in the database for various research programs. Our 
research protocol was approved by the Human Ethics Committee 
of the Japanese Red Cross Kumamoto Health Care Center 
(registration number: 137). 

Study population 

Healthy nonsmokers in the general population were sorted from 
those undergoing a medical checkup at the Japanese Red Cross 
Kumamoto Health Care Center during the two years from April 
2008 to March 2010. The medical checkup included a question- 
naire, test of physical strength and fitness, spirometry, chest X-ray, 
electrocardiogram (ECG), various blood tests, and medical 
examination by a physician. Nonsmokers were defined as the 
individuals who declared, in the questionnaire, that they never 
smoked before the medical checkup. Among the nonsmokers, the 
subjects who were reported to have no occupational history 
exposed to either biomass fuels or dusts and no obvious respiratory 



symptom, including dyspnea on exertion, nocturnal dyspnea, 
cough, sputum, or wheezing, were selected as the healthy 
nonsmoker candidates. The medical histories, as well as the results 
of various blood examinations, chest X-ray, and ECG, of these 
candidates were carefully inspected by the staff physicians. The 
candidates, who were confirmed to have neither cardiovascular 
disease nor respiratory diseases, such as lung cancer, bronchial 
asthma, COPD, interstitial lung disease, or infiltrative lung 
disease, were defined as healthy nonsmokers. The staff physicians 
excluded the nonsmokers who took the special medicines for 
asthma, COPD, and cardiovascular diseases excepting hyperten- 
sion. In addition, they excluded the nonsmokers having serious 
systemic diseases such as malignancy in any organ, renal failure 
requiring dialysis, diabetes mellitus with insulin therapy, and so on. 
However, the nonsmokers with hypertension or hyperlipidemia 
were accepted as the healthy nonsmokers unless their diseases were 
serious. Subsequently, the respiratory specialists checked the 
acceptability of spirometric data and finally decided the partici- 
pants eligible for the analysis. Through this process, the healthy 
nonsmokers with no spirometric data or incomplete and/or 
unacceptable traces of forced expiratory flow-volume curves were 
excluded. Thus, the total of 56,829 subjects over 20 years old 
(men: 34,160, women: 22,669) undergoing the medical checkup 
during the two years were assessed for eligibility and the 16,919 
healthy nonsmoking adults (men: 6,116, women: 10,803) were 
finally enrolled for the analysis. 

These participants were partitioned into three groups stratified 
by age at intervals of 20 years. The participants whose age was less 
than 39 years old were categorized as the young-age group, 
whereas those ranging between 40 to 59 years old were 
categorized as the middle-age group. The participants whose age 
was more than 60 years old were assigned to the advanced-age 
group. Among the men, the young-age, middle-age, and 
advanced-age groups consisted of 1,106, 3,343, and 1,667 subjects, 
respectively (Table 1). The number of female subjects forming the 
young-age, middle-age, and advanced-age groups was 1,578, 
6,427, and 2,798, respectively (Table 2). These classifications were 
made on the ground of the preliminary results obtained for the 12 
groups stratified by gender and age at intervals of 10 years, i.e. the 
group with age ranging from 20 to 29 years, 30 to 39 years, 40 to 
49 years, 50 to 59 years, 60 to 69 years, or more than 70 years in 
each gender. We estimated the effects of explanatory variables, 
including age and anthropometric factors of height, body weight, 
and body fat mass, on various spirometric parameters in these 12 
groups. The results showed that, in each gender, the effects of 
explanatory variables on spirometric parameters were qualitatively 
and quantitatively the same between the two successive age-groups 
in a first approximation (data not shown), indicating that the two 
successive age-groups would be united. Thus, the combination of 
the age-group with 20-29 years and that with 30-39 years formed 
the young-age group, whereas the group with 40-49 years and 
that with 50-59 years formed the middle-age group. Similarly, the 
unification of the group with 60—69 years and that with more than 
70 years made up the advanced-age group. We believe that the 
simplification of the group classification may hasten the compre- 
hension on the intricate results observed in this study. 

Spirometry and measurements of anthropometric 
measurements 

Forced expiratory pulmonary function tests were performed 
using an electric spirometer (DISCOM-21 FX, CHEST Co., 
Tokyo, Japan). Maneuvers were performed according to the 
standardization of lung function testing recommended by the 
ATS/ERS Task Force [20]. Under supervision by a skilled 
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Table 1. Demographic, anthropometric, and spirometric characteristics of male participants. 





Young-age Group (n = 1,106) 


Middle-age Group (n = 3,343) 


Advanced-age Group (n = 1,667) 


Age (years) 


35.7±2.8 


49.4±5.8" 


66.7±5.7*' t 


Height (cm) 


171.6+5.7 


169.8±6.0" 


164.7±5.9*' t 


Body weight (kg) 


70.0+11.1 


69.6±10.5 


63.2±8.5*' t 


BMI (kg/m 2 ) 


23.8±3.5 


24.1 ±3.2' 


23.3 ± 2.7*' 1 ' 


%FAT (%) 


22.9 ± 5.4 


22.8 ± 5.1 


20.5 ± 4.7*+ 


BFM (kg) 


16.5 ± 6.3 


16.2 ± 5.9 


13.2 ± 4.4* + 


FVC (L) 


4.25 ± 0.79 


3.95 ± 0.74* 


3.56 ± 0.67* 1 ' 


FEV, (L) 


3.51 ± 0.68 


3.17 ± 0.60* 


2.79 ± 0.55* 1 ' 


FEV,/FVC (%) 


82.72 ± 5.75 


80.50 ± 5.48* 


78.56 ± 5.67*' f 


PEF (L/s) 


8.41 ± 1.77 


8.19 ± 1.87* 


7.53 ± 1 .71 


FEF 50 (L/s) 


4.44 ± 1.34 


4.00 ± 1.27* 


3.47 ± 1 .21 


FEF 75 (L/s) 


1.58 ± 0.66 


1.20 ± 0.50* 


0.90 ± 0.44* 1 ' 



Values are the means ± SD. BMI: body mass index. %FAT: fat percentage of body mass. BFM: body fat mass calculated from body mass (kg) multiplied by %FAT/100. *: 
significantly different from the young-age group (at least p<0.01). + : different from the middle-age group (at least p<0.01). 
doi:1 0.1 371 /journal.pone.01 00733.t001 



technician for lung function tests, the participant repeated the 
forced expiratory maneuver until at least three acceptable and 
reproducible blows were obtained. Checking the values of FVC, 
FEVi, and extrapolated volume as well as the shape of the flow- 
volume curve, the technician selected the best FVC blow among 
the three traces. The trace thus selected was inspected by a 
respiratory specialist, who finally decided whether the best FVC 
trace of the participant was acceptable for the further analysis. We 
adopted FVC, FEVi, FEV, /FVC, PEF, FEF 50 , and FEF 75 as 
forced expiratory spirometric parameters for the subsequent 
analysis. 

A variety of anthropometric variables, including standing height 
(H, cm), body weight (BW, kg), body mass index (BMI), and fat 
percentage of body mass (%FAT, %), were examined for all 
participants. Standing height was measured with a stadiometer 
without shoes, while the participant was weighed on a scale 
wearing a suite of light clothes (weight: 500 g). Blood examinations 



and %FAT measurements were performed under fasting condi- 
tions without taking the breakfast. %FAT was measured with a 
bioelectrical impedance method (BF-220, TANITA Co., Tokyo, 
Japan). Body fat mass (BFM, kg) was calculated by multiplying BW 
by (% FAT/ 100). 

Multivariate analysis and selection of anthropometric 
variables 

Extending the model proposed by the GLI [3-5], we assumed 
that the effect of various explanatory variables (Ei, E 2 , E 3 ---E n ) on 
the spirometric parameter (SP) was described by the multiplicative 
exponential function. Thus, the regression model for a certain SP 
assumed the following: 

SP = K 0 -(Ei) al -(E 2 ) a2 -(E 3 ) a3 -(E n ) an (1) 
where Kn is the constant equal to Exp(a 0 ). By taking the logarithm 



Table 2. Demographic, anthropometric, and spirometric characteristics of female participants. 



Young-age Group (n = 1,578) Middle-age Group (n = 6,427) Advanced-age Group (n = 2,798) 



Age (years) 


35.6 ± 3.0 


49.7 ± 5.8* 


66.1 ± 5.4*' 1 ' 


Height (cm) 


159.5 ± 5.2 


156.9 ± 5.4* 


152.3 ± 5. I* 1 ' 


Body weight (kg) 


53.3 ± 8.6 


54.4 ± 8.6* 


52.1 ± 7.6*' 1 ' 


BMI (kg/m 2 ) 


20.9 ± 3.2 


22.1 ± 3.3* 


22.4 ± 3.0*' 1 ' 


%FAT (%) 


26.2 ± 6.4 


27.8 ± 6.2* 


27.9 ± 5.8* 


BFM (kg) 


14.4 ± 6.1 


15.6 ± 6.0* 


14.9 ± 5.1* f 


FVC (L) 


3.25 ± 0.59 


3.12 ± 0.64* 


2.82 ± 0.72* 1 ' 


FEV, (L) 


2.73 ± 0.46 


2.53 ± 0.51* 


2.23 ± 0.59* 1 ' 


FEV,/FVC (%) 


84.27 ± 6.32 


81.28 ± 5.65* 


79.24 ± 5.47*' 1 ' 


PEF (L/s) 


6.01 ± 1.49 


6.09 ± 1.56 


5.65 ± 1 .65* f 


FEF 50 (L/s) 


3.60 ± 0.96 


3.31 ± 1.01* 


2.87 ± 1.08* f 


FEF 75 (L/s) 


1.35 ± 0.55 


0.99 ± 0.44* 


0.75 ± 0.42* 1 ' 



Values are the means ± SD. Symbols are the same as denoted in Table 1 . *: different from the young-age group (at least p<0.03). ': different from the middle-age group 
(at least p<0.03). 

doi:1 0.1 371 /journal.pone.01 00733.t002 
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(Ln) of both sides of eq. (1), the model is converted to the 
logarithmic additive function as follows: 

Ln(SP) = a 0 + ai -LnfE, ) + a 2 -Ln(E 2 ) 

(2) 

+ a 3 -Ln(E 3 ) + --- + a n -Ln(E n ) 

In eq. (2), a ; (i = 1 to n) is the partial regression coefficient of a 
particular explanatory variable, whereas a 0 is the invariable 
constant. The most important issue is what kinds of anthropo- 
metric factors should be introduced as explanatory variables 
deciding SPs. In addition to age (A), height (H) has universally 
been accepted as the principal explanatory variable predicting 
SPs. Other anthropometric variables, such as body weight (BW), 
body mass index (BMI), body fat mass (BFM, defined as 
multiplying BW by fat percentage of body mass (%FAT)), fat- 
free mass (FFM), and body surface area (BSA), have also been 
taken as explanatory variables deciding SPs [6-19]. BW is 
important because it is the sum of the various constituents forming 
the body, including respiratory muscles that influence a variety of 
pulmonary function parameters [7-9,21]. BFM may also play a 
significant role because fat accumulation along central and 
peripheral airways elicits airway narrowing and leads to a decrease 
in various SPs [7,22]. FFM is important because respiratory 
muscles are contained in FFM. In the present study, however, we 
did not measure FFM. This is because, for estimating FFM, the 
skinfold thickness at various regions of the body should be 
measured [23], the method being time-consuming and not fitted 
for medical checkups subject to many persons. Instead, we 
assumed that the difference between BW and BFM might act as 
the component reflecting FFM in a first approximation. Although 
some authors used BMI and BSA as explanatory variables for 
predicting SPs [7-9,10,11,18], we considered that the effects of 
these two variables on SPs would be involved in those of H and 
BW. This is due to the fact that BMI and BSA are expressed as 
BW/H 2 and c 0 -(BW) <:1 -(H) <:2 (c ( „ c b and c 2 : constants), respec- 
tively. Taking the logarithm of BMI or BSA, each value is 
converted to the logarithmic additive function as: 

Ln(BMI) = Ln(BW)-2-Ln(H) (3) 



Ln(BSA) = Ln(c 0 ) + C! -Ln(BW) + c 2 -Ln(H) (4) 

Eqs. (3) and (4) certainly indicate that the effect of BMI or BSA 
on SPs is partitioned into the effects of BW and H. Based on these 
reasons, we adopted age (A), height (H), body weight (BW), and 
body fat mass (BFM) as explanatory variables for predicting 
reference means of various spirometric parameters in the present 
study. Thus, the regression model for a certain SP was assumed 
the following: 

Ln(SP) = a 0 + ai -Ln(A) + a 2 -Ln(H) 
+ a 3 -Ln(BW) + a 4 -Ln(BFM) 

We decided the coefficients of a; (i = 0 to 4) in each age group of 
either gender by applying the classical multiple-regression analysis 
with least-squares minimization. Overall agreement between the 



predicted reference means and observed values was judged by the 
coefficient of determination adjusted for degrees of freedom 
(adjusted-R 2 ). Normality for the distribution of residuals of each 
Ln(SP) was examined using the Kolmogorov-Smirnov and 
Shapiro-Wilk tests. These tests indicated that the residuals between 
the logarithm-transformed forms of observed values and those of 
the predicted reference means did not follow the normal 
distribution for most of the spirometric parameters. Therefore, 
we did not calculate the lower limit of normal (LLN) of each 
Ln(SP) from the standard deviation of residuals (RSD). Instead, we 
determined it from the 5 th percentile of distribution of the 
residuals. The disparity between the log-transformed value of 
predicted reference mean and LLN was defined as ALLN as 
follows: 

ALLN = log-transformed reference mean — LLN (6) 

All calculations were performed using the IBM SPSS Statistics 
(Version 21.0, SPSS Inc., an IBM Co., NY, USA). Unless 
otherwise specified, the values were expressed as the means ± 
standard deviations. A p-value lower than 0.05 was deemed to be 
statistically significant. 

Results 

In the male participants, H was the maximum in young-age 
group, whereas BMI displayed the largest value in the middle-age 
group (Table 1). The values of BW, %FAT, and BFM in young- 
age and middle-age groups were higher than those in the 
advanced-age group. In the men, all spirometric parameters were 
significantly different between three groups, and the highest value 
was observed in the young-age group (Table 1). 

Somewhat different from the male participants, H was 
decreased and BMI increased with aging in the female participants 
(Table 2). BW or BFM reached the maximum level in the middle- 
age group, whereas %FAT displayed the highest value in the 
middle-age and advanced-age groups. In the female participants, 
spirometric parameters other than PEF revealed trends identical to 
those investigated in the male participants (Table 2). 

Among the partial regression coefficients of a 0 , a 1; a 2 , a 3 , and a 4 
determined in the young-age, middle-age, and advanced-age 
groups of both genders, only the values reaching statistical 
significance are depicted in Tables 3, 4, and 5, respectively. 

Relationship between the anthropometric variables and 
spirometric parameters of the young-age group (Table 3) 

Age did not act as a significant determinant on FVC and PEF in 
both genders as well as on FEV! in male participants. However, 
age was a decisive factor for FEVi in female participants and 
FEV] /FVC as well as flow parameters, including FEF 50 and 
FEF75, in both genders. 

H had a substantial impact on FVC, FEV], and the flow 
parameters of FEF 50 and FEF 75 in both genders. However, it did 
not have a significant impact on FEVi/FVC in both genders and 
PEF in female participants. 

BW functioned as one of the determinants on FEV! /FVC in 
male participants and FVC, FEVi PEF, and FEF 50 in female 
participants. 

BFM had an appreciable impact on FEV! /FVC in female 
participants and FEF 73 in both genders. However, there was no 
spirometric parameter simultaneously influenced by BW and 
BFM. 
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Figure 1. Age-specific effects of height (H) on decision of FVC 
and FEV, estimated for three different age groups of either 
males or females. (A): Partial regression coefficients of Ln(H) for 
reference means of Ln(FVC) in young-, middle-, and advanced-age 
groups of both genders. Partial regression coefficients of Ln(H) for 
Ln(FVC) are denoted as ALn(FVC)/ALn(H). (B): Partial regression 
coefficients of Ln(H) for Ln(FEV q ) are designated as ALn(FEV,)/ALn(H). 
doi:10.1371/journal.pone.0100733.g001 
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Figure 2. Age-specific impacts of height (H) on PEF and FEF 50 
estimated for three different age groups of both men and 
women. (A): Partial regression coefficients of Ln(H) for reference means 
of Ln(PEF) in young-, middle-, and advanced-age groups of either 
gender. Partial regression coefficients of Ln(H) for ALn(PEF) are defined 
as ALn(PEF)/ALn(H). (B): Partial regression coefficients of Ln(H) for 
Ln(FEF 50 ) are defined as ALn(FEF 50 )/ALn(H). 
doi:1 0.1 371/journal.pone.01 00733.g002 



Middle-age group (Table 4) 

Unlike the young-age group, age played an important role in 
deciding all spirometric parameters in the middle-age group 
irrespective of gender. 

Although H acted as a significant explanatory variable 
determining all spirometric parameters in male participants, it 
had no significant influence on the flow parameters of PEF, FEF 50 , 
and FEF 75 in female participants. 

BW and BFM concurrently exerted a considerable influence on 
FEVi and FEF 75 in male participants as well as FEVi/FVC and 
all flow parameters in female participants. The direction of 
contributions of BW and BFM was opposite, i.e., BW increased, 
whereas BFM decreased, a given spirometric parameter. However, 
apart from BW, BFM appreciably contributed to determination of 
the FVC and FEVi /FVC of male participants. 

Advanced-age group (Table 5) 

Similar to the middle-age group, age played a significant role for 
all spirometric parameters independent of gender in the advanced- 
age group. It should be noted that the decline of pulmonary 
function with age was consistently larger in the advanced-age 
group than in other age groups regardless of the spirometric 
parameter or the gender. 

H functioned as a decisive factor for many spirometric 
parameters in advanced-age group, but it had no significant 



influence on FEVi /FVC m male participants and FEF 50 in female 
participants. 

In contrast to the middle-age group, the concurrent contribu- 
tion of BW and BFM to the spirometric parameters was not 
investigated in the advanced-age group. BW decreased FEF 75 in 
male participants but increased FEVi /FVC an d FEF50 in female 
participants. BFM increased FVC in male participants as well as 
PEF in female participants but decreased FEVi m male 
participants. 

Age- and gender-dependent contributions of 
anthropometric variables to spirometric parameters 
(Figs. 1 to 9) 

Depending on gender and age, the anthropometric variables of 
H, BW, and BFM were inextricably linked with spirometric 
parameters in a qualitatively and/or quantitatively different 
manner. One exception was the effect of H on FVC and FEVi. 
For these parameters, the impact of H was relatively homogeneous 
and approximately independent of age in both genders (Fig. 1). 
The identical tendency, i.e., the constant partial regression 
coefficient of H with respect to age, approximately held for PEF 
and FEF50 but not for FEV, /FVC and FEF 75 for the male 
participants (Figs. 2, 3). On the other hand, the partial regression 
coefficients of H for spirometric parameters, except for FVC and 
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Figure 3. Age-specific effects of height (H) on decision of FEV,/ 
FVC and FEF 75 estimated for three different age groups of 
male and female participants. (A): Partial regression coefficients of 
Ln(H) for reference means of Ln(FEV q /FVC) in young-, middle-, and 
advanced-age groups of either gender. Partial regression coefficients of 
Ln(H) for Ln(FEV,/FVC) are denoted as ALn(FEV,/FVC)/ALn(H). (B): Partial 
regression coefficients of Ln(H) for Ln(FEF 75 ) are designated as 
ALn(FEV 75 )/ALn(H). 

doi:10.1371/journal.pone.0100733.g003 



FEV^ were not constant and changed significandy with age in the 
female participants (Figs. 1-3). 

BW exerted little influence on FVC, PEF, and FEF 30 in any of 
the age groups in the male participants (Figs. 4, 5). For the male 
participants, however, BW had appreciable impacts on FEVi, 
FEV! /FVC and FEF 75 in some age groups (Figs. 4, 6). In the 
female participants, BW played an important role in determining 
FVC or FEVi only in the young-age group (Fig. 4). The FEV \/ 
FVC and the flow parameters of PEF, FEF 50 , and FEF 75 in the 
female groups were substantially influenced by BW in an age- 
dependent manner, with the maximum effect being detected in the 
middle-age group (Figs. 5, 6). 

Although H and BW generally were factors that increased 
spirometric parameters, BFM exerted negative impacts on most of 
the spirometric parameters in an age-dependent fashion (Figs. 7, 
8). The concurrent impacts of BW and BFM were specifically 
detected in the female middle-age group, in which BW functioned 
as an increasing factor, whereas BFM was as a decreasing factor 
for a variety of flow parameters (Fig. 9). The influence of BFM on 
spirometric parameters was weaker in the young- and advanced- 
age groups irrespective of gender (Fig. 7, 8). 
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Figure 4. Age-specific contributions of body weight (BW) to 
FVC and FEV 1 estimated for three different age groups of both 
genders. (A): Partial regression coefficients of Ln(BW) for reference 
means of Ln(FVC) in young-, middle-, and advanced-age groups of 
either gender. They are designated as ALn(FVC)/ALn(BW). (B): Partial 
regression coefficients of Ln(BW) for Ln(FEV,) are designated as 
ALn(FEV,)/ALn(BW). 
doi:10.1 371/journal.pone.01 00733.g004 

Discussion 

Critique of the method 

Although we applied the multiplicative model similar to that 
proposed by the GLI [3-5], we did not determine the partial 
regression coefficients for Ln(SP) in terms of the LMS method with 
smoothing such as penalized beta-spline [24,25], i.e., more 
generalized additive, linear model describing location, scale, and 
shape (GAMLSS) of a given dataset. Instead, we determined the 
coefficients with the method of classical multiple-regression 
analysis of least-squares minimization without smoothing. The 
reason for this simplification is that the main purpose of the 
present study is not to establish the generalized regression 
equations for SPs applicable for all ages but to explore the age- 
and gender-associated impacts of each anthropometric factor on 
each SP. 

We should note that it has not been verified whether the age- 
related importance of anthropometric variables for determining 
spirometric parameters investigated in the present study holds true 
for all races and ethnicities because our study population consisted 
entirely of nonsmoking healthy Japanese adults. 

Aging- and gender-dependent contributions of age to 
spirometric parameters 

Although age was not a determinant on FVC and PEF in the 
young-age groups of both genders and FEV] in the male young-age 
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Figure 5. Age-specific impacts of body weight (BW) on PEF and 
FEF 50 estimated for three different age groups of male and 
female participants. (A): Partial regression coefficients of Ln(BW) for 
the reference means of Ln(PEF) in young-, middle-, and advanced-age 
groups of both genders. Partial regression coefficients of Ln(BW) for 
Ln(PEF) are defined as ALn(PEF)/AI_n(BW). (B): Partial regression 
coefficients of Ln(BW) for Ln(FEF 50 ) are defined as ALn(FEF 50 )/ALn(BW). 
doi:10.1371/journal.pone.0100733.g005 

group, age played an important role in deciding other spirometric 
parameters in any age groups irrespective of the gender. It should be 
noted that the decline of pulmonary function with age was 
consistently larger in the advanced-age group than in other age 
groups regardless of the spirometric parameter or the gender, 
indicating that the degree of effects of age on spirometric parameters 
would differ depending on aging. These facts suggest that the 
assumption of aging-independent, constant partial regression 
coefficients of age are not acceptable for any spirometric parameters 
in both genders. 

Age- and gender-dependent contributions of 
anthropometric variables to spirometric parameters 

We demonstrated that the anthropometric variables of H, BW, 
and BFM were linked with spirometric parameters generally in an 
age-dependent manner. However, the impact of H on FVC and 
FEVi was relatively homogeneous and approximately indepen- 
dent of age in both genders, indicating that the assumption of an 
age-independent, constant value of the partial regression coeffi- 
cient of H is only acceptable for the equation predicting the 
reference mean of FVC or FEVi irrespective of the gender. The 
identical assumption, i.e., the constant partial regression coefficient 
of H with regard to age, approximately held for PEF and FEF 50 
for the male participants. On the other hand, the partial regression 
coefficients of H for most of the spirometric parameters in the 
female participants were changed significantly with age, suggesting 
that the age-dependent contribution of H to the spirometric 
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Figure 6. Age-specific contributions of body weight (BW) to 
decision of FEV,/FVC and FEF 75 estimated for three different 
age groups of both genders. (A): Partial regression coefficients of 
Ln(BW) for reference means of Ln(FEV,/FVC) in young-, middle-, and 
advanced-age groups of both genders. They are designated as 
ALn(FEV,/FVC)/ALn(BW). (B): Partial regression coefficients of Ln(BW) 
for Ln(FEF 75 ) are designated as ALn(FEF 75 )/ALn(BW). 
doi:10.1 371/journal.pone.01 00733.g006 

parameters is more evident in the female participants than in the 
male participants. 

For the male participants, BW exerted little influence on FVC, 
PEF, and FEF50 in any of the age groups, allowing us to set the 
partial regression coefficients of BW for these spirometric 
parameters to zero, i.e., the homogeneous effects of BW on 
FVC, PEF, and FEF 50 in the male participants. For the male 
participants, however, BW had appreciable impacts on FEVi, 
FEVi /FVC an d FEF75 in some age groups, indicating that the 
assumption of the age-independent, constant partial regression 
coefficients of BW was not acceptable for these spirometric 
parameters. For the female participants, BW played a substantial 
role in determining all of the spirometric parameters in an age- 
dependent manner. Joining the findings observed for male and 
female participants, we considered that, although the impact of 
BW on spirometric parameters were very inhomogeneous and 
changed distincdy depending on gender and age, they were more 
conspicuous for the female participants than for the male 
participants. 

In contrast with H and BW, BFM exerted negative influences 
on most of the spirometric parameters in an age-dependent 
fashion. In the female middle-age group, BW functioned as an 
increasing factor, whereas BFM was as a decreasing factor for a 
variety of flow parameters, indicating that the negative effect of 
BFM on certain spirometric parameters could be offset by the 
positive effect of BW. The importance of BW and BFM in the 
female middle-age group is supported by the fact that physique 
changes with aging, leading to both BW and BFM reaching the 
maximums in this female age group (Table 2). On the other hand, 
the influence of BW and BFM on spirometric parameters 
appeared to be weak in the advanced-age population irrespective 
of gender. This is also explained by the finding that BW and BFM 
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Figure 7. Age-specific impacts of body fat mass (BFM) on FEV, 
and FEV-i/FVC estimated for three different age groups of both 
men and women. (A): Partial regression coefficients of Ln(BFM) for 
reference means of Ln(FEV q ) in young-, middle-, and advanced-age 
groups of either gender, which were expressed as ALn(FEV,)/ALn(BFIVI). 
(B): Partial regression coefficients of Ln(BFM) for Ln(FEV,/FVC), being 
defined as AI_n(FEV q /FVC)/ALn(BFIVI). 
doi:10.1371/journal.pone.0100733.g007 

generally decreased with aging, reaching the minimums or the 
lower levels in the advanced-age population (Tables 1, 2). 

Importance of body weight and body fat mass in 
determining spirometric parameters 

BW is the sum of the various constituents, including respiratory 
muscles and fat deposition in the thoracic cavity and around the 
airway. Respiratory muscles influence the maximal respiratory 
pressures, and, hence, a variety of pulmonary function parameters 
such as inspiratory capacity, FVC, FEVi, and PEF [7-9,21]. Fat 
components represented by BFM influence the lung volumes, the 
work of breathing, and, in some circumstances, the airway caliber 
[7,22]. Fat accumulation along a large airway may cause central 
airway narrowing, leading to a decrease in FVC, FEVi, and/or 
PEF. On the other hand, fat deposition along a small airway may 
elicit peripheral airway narrowing, leading to reduction in flow 
parameters such as FEF 50 and FEF 75 . We found that the negative 
impact of BFM on the spirometric parameter was age-dependent 
and evident for FVC (male), FEV[ (male), FEV^FVC (both 
genders), PEF (female), FEFr, 0 (female), or FEF 75 (both genders) in 
the middle-age group. Our findings are partially consistent with 
those reported by Cotes et al. [9], who demonstrated that the fat 
percentage of body mass (%FAT) plays a negative role in 
determining FVC and FEV! in both genders. Unfortunately, 
however, they did not address the age-dependent contribution of 
%FAT to the spirometric parameters. Interestingly, we found that, 
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Figure 8. Age-specific impacts of body fat mass (BFM) on 
deciding FEF 50 and FEF 75 estimated for three different age 
groups of male and female participants. (A): Partial regression 
coefficients of Ln(BFM) for reference means of Ln(FEF 50 ) in young-, 
middle-, and advanced-age groups of either gender. They are expressed 
as ALn(FEF 50 )/ALn(BFM). (B): Partial regression coefficients of Ln(BFM) 
for Ln(FEF 75 ) are defined as ALn(FEF 75 )/ALn(BFM). 
doi:1 0.1 371/journal.pone.01 00733.g008 



in contrast to BFM, BW had a positive impact on most of the 
spirometric parameters, though BFM was one of the components 
forming BW. This peculiar finding may be explained, in part, by 
the fact that the content of respiratory muscles is among the 
components of BW as well. Our findings are in accordance with 
the study of Pistelli et al. [18], who demonstrated the positive effect 
of BMI on many spirometric parameters, but our findings are 
seemingly inconsistent with that of Chen et al. [6] and Chinn et al. 
[7], both of whom demonstrated a weight-gain-associated decline 
in lung function. A longitudinal study done by Chen et al. [6] 
showed that each kg of weight-gain was associated with an excess 
loss of 26 mL in FVC and 23 mL in FEVi in men, and 14 mL and 
9 mL respectively in women. Chinn et al. [7] revealed that FVC 
and FEVi were decreased by about 20 mL per kg-gain of BW in a 
longitudinal study for male shipyard workers. As pointed out by 
Chinn et al. [7] , however, the increase in BW of a certain person 
in a longitudinal observation is generally caused by the increase in 
BFM but not by the increase in FFM. In other word, the effect of 
weight-gain on lung function in a longitudinal observation reflects 
the effect of increased BFM at a relatively constant level of FFM. 
Therefore, we considered that the findings reported by Chen et al. 
[6] and Chinn et al. [7] were not inconsistent with our findings, in 
which BFM was found to have a negative impact on lung function. 
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Figure 9. Opposite effects of body weight (BW) and body fat mass (BFM) on FEF 75 in female middle-age group. (A): Ln(FEF 75 ) vs. 
Ln(BW). Parentheses: absolute values of BW without log-transformation. Red line: slope of the equation predicting reference means of Ln(FEF 75 ) from 
Ln(BW), which is equal to partial regression coefficient shown in Table 4 and Fig. 6-(B). (B): Ln(FEF 75 ) vs. Ln(BFM). Parentheses: absolute values of BFM 
with no log-transformation. Red line: slope of the equation predicting reference means of Ln(FEF 75 ) from Ln(BFM), which is equal to partial regression 
coefficient given in Table 4 and Fig. 8-(B). 
doi:1 0.1 371 /journal. pone.0100733.g009 



Clinical implication of the present Study means of various spirometric parameters. Furthermore, we 

We found that, in addition to age and height, body weight and demonstrated that the effects of all explanatory variables introduced 
body fat mass had substantial impacts on deciding the reference in the present study were conspicuously age-dependent. Therefore, 
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we considered that our findings would serve as a warning to the 
reference equations reported over several decades, in which aging- 
independent contributions of two explanatory variables, including 
age and height, to spirometric parameters were simply assumed. 
Since spirometric parameters are measured under a forced 
expiratory maneuver, the effect of the variables involving the 
respiratory muscles should be taken into account. For this purpose, 
fat-free mass (FFM) or bone-free lean body mass (BF-LBM) deserves 
a good indicator [17]. It is not easy, however, to measure these 
variables in a cohort consisting of numerous subjects. Therefore, we 
used body weight as a substitute for them, in which respiratory 
muscles are contained. However, we should notice that body weight 
has a disadvantage, i.e., it includes body fat mass (BFM) as well. In 
an obese subject, the excessive body weight caused by increased 
BFM will conceal the positive effect of respiratory muscles in 
association with disclosing the negative effect of fat deposition in the 
thoracic and abdominal cavities, which leads to narrowing the 
airway calibers and inhibiting the diaphragm movement. We 
consider that, to defeat the disadvantageous effect of body weight 
caused by increased BFM, it is indispensable to concurrently 
estimate the independent effects of body weight and that of BFM, 
like the present study did. 

Although we did not estimate the quantitative difference 
between the predicted means of spirometric parameters calculated 
from our equations (aging-dependent differences in partial 
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Conclusions 
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